Thermoelectric detection of ferromagnetic resonance of a nanoscale ferromagnet 
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We present thermoelectric measurements of the heat dissipated due to ferromagnetic resonance of 
a Permalloy strip. A microwave magnetic field, produced by an on-chip coplanar strip waveguide, 
is used to drive the magnetization precession. The generated heat is detected via Seebeck mea- 
surements on a thermocouple connected to the ferromagnet. The observed resonance peak shape is 
in agreement with the Landau-Lifshitz-Gilbert equation and is compared with thermoelectric finite 
element modeling. Unlike other methods, this technique is not restricted to electrically conductive 
media and is therefore also applicable to for instance ferromagnetic insulators. 
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Thermal effects in ferromagnetic materials are sub- 
ject to extensive research since the discovery of the 
spin-Seebeck effect [THS]. Recently, spin dynamics and 
(spin-) caloritronics, two popular branches of spintron- 
ics, started to come together as spin pumping induced 
by spin dynamics has been proposed as the origin of the 
spin-Seebeck effect [H |3] . Magnetization dynamics has 
been studied thoroughly in magnetic systems as it is an 
important mechanism for future spintronic applications, 
e.g. for microwave generators [SI [7] and spin sources via 
spin pumping [51 [5]. However, dissipation mechanisms 
that accompany magnetization dynamics, and cause lo- 
cal heating, are still not fully understood [TOl [H] . 

Here, we focus on a new aspect, the coupling between 
magnetization dynamics and the generation of heat. We 
deduce from thermoelectric measurements on a Permal- 
loy (Py) island the heat dissipation during ferromagnetic 
resonance. This on-chip detection technique, based on 
the Seebeck effect, offers a novel method for characteriz- 
ing ferromagnetic resonance and hence is distinctly differ- 
ent from other techniques, such as scanning thermal mi- 
croscopy [T21[T3|. Due to the thermal detection, electrical 
contact to the ferromagnet is in principle not required. 
Hence, this method allows for FMR measurements on 
non-conductive materials like ferromagnetic insulators. 

When a ferromagnet is brought into resonance, energy 
is absorbed from the applied microwave field. This en- 
ergy causes the magnetization M to precess around an 
effective field H and the motion is well described by 
the Landau-Lifshitz-Gilbert (LLG) equation dM /dt — 
-jM xH+{a/Ms)MxdM/dt with = 176 GHz/T the 
gyromagnetic ratio. The last term in the LLG equation 
describes the damping of the magnetization towards the 
direction of the effective field H, using the phenomeno- 
logical damping parameter a. This process is purely 
dissipative and converts magnetostatic energy into heat. 
During ferromagnetic resonance (FMR) this continuous 
dissipation leads to heating of the ferromagnetic mate- 
rial. 

In this experiment, we measured the temperature of a 
ferromagnet while subject to a microwave magnetic field. 



A ferromagnetic strip is placed close to the shortened 
end of a coplanar strip waveguide (CSW) as shown in 
Fig. [l] Microwave power is applied to the CSW, leading 
to an out of plane rf magnetic field. A static magnetic 
field is applied along the easy axis of the magnet. In 
addition, a thermocouple consisting of a NiCu and Pt 
wire is connected to the ferromagnet by a Au bridge. In 
this way, the temperature can be measured by making 
use of the Seebeck effect. The Seebeck effect describes 
the generation of a voltage due to a temperature gradi- 
ent, Vy = —SWT, with S the material dependent See- 
beck coefficient. The voltage that develops across the 
Pt wire is different than the voltage that develops across 
the NiCu wire, leading to a nonzero voltage between the 
two wires. This thermovoltage scales with the Seebeck 
coefficients (S'pt — S'NiCu) and the temperature difference 
(Ti — To). Note that NiCu is chosen because of its rel- 
atively high Seebeck coefficient (S'NiCu — -'iO^iV/K and 
Spt = -5^V/K). In the following, we calculate the tem- 




FIG. 1. (color online) Concept of the thermoelectric detec- 
tion of ferromagnetic resonance. A coplanar strip waveguide 
generates a microwave magnetic field h^f in the z direction 
acting on a small ferromagnetic strip. A static magnetic field 
ho is applied along the x axis. When the strip is brought into 
resonance, it absorbs energy from the field which is dissipated 
as heat. The dissipation is detected by a Pt-NiCu thermocou- 
ple which probes the temperature Ti of the ferromagnet with 
respect to a reference temperature To via the Seebeck effect. 
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perature rise during FMR from the dissipated power. 

The energy of a ferromagnetic particle in a magnetic 
field, the Zeeman energy, is given by: 



E = - M ■ BdV 
'v 



(1) 



with M the magnetization, B = fio (^-ffcxt + Hd /2 j the 
sum of the externally applied magnetic field and the de- 
magnetizing field and V the volume of the particle ' 14j . 
Here, H]j is divided by two to compensate for counting 
each volume twice in the integral. For this experiment, a 
static magnetic field, /ig, is applied in the x direction and 
a driving rf field, hrf cos cut, in the z direction, making 
B = {ho—Na:mx/2,—Nymy/2,hrfCosujt — Nzmz/2) 
with Nx, Ny and Nz the demagnetization factors. The 
dissipation energy can now be calculated from the time 
derivative of i?, assuming a uniform M and B: 



dE / dM ^ dB 

dt \ dt dt 



(2) 



where the first part of Eq. [2]expresses the dissipation due 
to the magnetization motion and the second part the en- 
ergy absorbed from the microwave field. In equilibrium, 
the absorption of energy from the microwave field equals 
the dissipation, (dE/dt) — 0, leading to heating of the 
ferromagnet. In order to find an expression for the dis- 
sipated power, we use a procedure similar to Ref. [15] 
where the magnetization dynamics is described by the 
linearized LLG equation. We assume that for small angle 
precessional motion dm^/dt = such that ttIx is constant 
and the solution to the Landau-Lifshitz-Gilbert (LLG) 
equation can be written in terms of the sum of in-phase 
and out-of-phase susceptibilities. The components my 
and niz are now defined as niy = x'yUJi cosut+Xy^i sinwt 
and niz = Xz^i cosujt + Xz'-^i simut with: 



auj'^{ujy + Wz) 



[uj'^ — LOyLOzY + {aU)Y[0Jy + LOzY 
II _ - i^yi^z) 

y (a;2 - ujyUJzY + {aujY{u:y + uOzY 



(3) 



Here, uj is the frequency of the driving rf field, ujy — 
7(/io - {N.^ - Ny)ms), LUz = -/{ho - {No, - Nz)ms), 
LJi — jhrf and the saturation magnetization. The 
x'z and x'z are related via x'z = '^x'y — ^y/^Xy and 
x'z = ctXy + '^y/'^Xyj respectively. With these expres- 
sions for rrix, rUy, rriz and B^, By, Bz one can easily find 
its time derivatives and calculate the relevant dot prod- 
uct of Eq. [2j We are not interested in high frequency 
variations in the dissipated power and hence, average out 
all contributions that vary with a frequency uj or 2uj and 
find: 



dE\ 

I dissipation 
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FIG. 2. (color online) a) Scanning electron microscope (SEM) 
image of the thermoelectric FMR device, b) Image of a device 
with four contacts for dc AMR detection of FMR. 



From this expression we can deduce that the resonance 
peak shape of the dissipated power is determined by 
x" and scales with the applied microwave field and fre- 
quency. In order to convert this power into a tempera- 
ture rise, we make use of 3D finite element thermoelectric 
modeling. For details about the modeling we refer to ear- 
lier publications [121 [T7] . 

The samples are fabricated using a three-step electron 
beam lithography process on top of a thermally oxidized 
Si substrate. A scanning electron microscope (SEM) im- 
age of the investigated devices is shown in Fig. [2] The 
devices consist of a 50nm thick Py strip (2 /im x 400 nm) 
close to a 100 nm thick Au coplanar strip waveguide. The 
CSW is made using an optical lithography process. In the 
thermoelectric device (Fig. [2^), there are two 40 nm thick 
contacts (Ft and NiCu) forming a thermocouple. The 
Py island is connected to the thermocouple by a highly 
thermal conductive Au contact (thickness: 120 nm). The 
other side of the thermocouple is connected by 120 nm 
thick Au contacts to the bonding pads. In the case of the 
anisotropic magnetoresistance (AMR) device (Fig. [2]d), 
four 120 nm thick Au contacts directly connect to the Py 
strip. The NiCu is deposited by DC sputtering to pre- 
serve the original alloy composition (45% Ni, 55% Cu). 
To avoid lift-off problems a double-layer resist technique 
with a large undercut (PMMA-MA and PMMA 950K) is 
used. The Au, Ft and Py are deposited using an e-beam 
evaporator (base pressure 1 x 10"'' mbar) and a single 
layer resist (PMMA 950K). Prior to the Au deposition, 
the NiCu, Ft and Py surfaces are cleaned with Ar ion 
milling. 

For the measurements, we have used a frequency mod- 
ulation method to obtain a better signal to noise ratio 
and to remove background voltages due to heating of 
the CSW short. The microwave field frequency is al- 
ternated between two different values with a separation 
of 5 GHz. A lock- in amplifier, tuned to the same fre- 
quency (17 Hz), measures the difference in dc voltage 
across contacts 1 and 2 (Fig. [2^) between the two fre- 
quencies iy = 14=high — Vi=\ovi). Because of the large 
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separation between T4=iow and Vf=high, they can not both 
fulfill the resonance condition at a specific magnetic field. 
With this method, one effectively measures the difference 
in the Seebeck or AMR voltage when the ferromagnet is 
in- and off-resonance. All measurements were performed 
at room temperature. 

Fig. [3^ shows the measured Seebeck voltage as a func- 
tion of magnetic field for different rf field frequencies (10 
- 20 GHz) for 12 dBm rf power. The position of the peaks 
and dips correspond to the resonance field for /low and 
/high, respectively. We have plotted the peak position 
as a function of the applied rf frequency in Fig. [3]3 and 
found peak heights ranging from 46 nV at 10 GHz to 
105 nV at 17 GHz. For a uniform precessional mode, the 
resonance field is related to w by the Kittel equation [TH] : 

= l^{ho - {N, - Ny)m,){hn ~ {N., - iV,)m,) (5) 

The line corresponds to a fit of Eq. [5] and the almost 
perfect fit confirms the uniform precessional motion. We 
obtained the following fitting parameters: = 0.01, 
Ny = 0.09, = 0.90 and fiorris = 1.11 T. These param- 
eters have been used to calculate the dissipated power of 
Eq. |4] for different frequencies. However, in order to do 
this accurately one first need to determine the magnitude 
of the rf magnetic field experimentally. 

To obtain the correct experimental value for the mag- 
nitude of the rf magnetic field, we have measured the dc 
anisotropic magnetoresistace (AMR) in a dedicated de- 
vice with a four-terminal geometry (shown in Fig. ^p) . 
The AMR effect describes the dependence of the resis- 
tance on the angle between the current / and the di- 
rection of magnetization M by i? = Ro—ARsin^ 9, where 
Rq is the resistance of the strip when / and M are par- 
allel, and AR the difference in resistance between the 
parallel and perpendicular alignment of / and M. A mea- 
surement of i? as a function of a perpendicular applied 
magnetic field is plotted in Fig. [4j:. From this measure- 
ment, we determined the magnitude of the AMR effect 
and found AR/Rq = 1.5 %. For a steady resonant pre- 
cession, the average cone angle 9c of the precession can 
now be extracted from the observed AMR voltage. 

Fig. displays the AMR voltage versus magnetic field 
for different microwave field frequencies. For this mea- 
surement, we have used a dc current I^c of 300 /iA. 
The obtained voltage now corresponds to the dc cur- 
rent multiplied with the resistance change, being V — 
/dcAi?sin^ 9c, and the precession angle can be extracted. 
Using Eq.[3]and an expression for the average cone angle: 

m-c.lix'y'+x'^' + X? + xT)/2 (6) 

one can deduce hrf = wi/7 from a fit of the measured 
peak height, and the result is plotted in Fig. |4f3 for a 
microwave power of 12 dBm. The field strength is found 
to decrease twofold when the frequency is increased from 
10 to 20 GHz. We attribute this to frequency dependent 
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FIG. 3. (color online) a) Series of Seebeck voltage versus 
magnetic field measurements for 11 difi^erent frequencies. The 
traces are offset by 150 nV for clarity. Due to the modula- 
tion technique using two driving frequencies that are 5 GHz 
apart, peaks and dips are observed at the resonance fields for 
both frequencies, b) Frequency versus the magnetic field at 
the center of the resonance peak. The line corresponds to 
a fit of the Kittel equation, c) Generated power and corre- 
sponding Seebeck voltage calculated using Eq. [4]and thermo- 
electric finite-element modeling for multiple frequencies. The 
measured peak heights of a) are indicated by the black dots. 



attenuation of the microwave signal, leading to smaller rf 
fields at higher frequencies. 

Now we can calculate, using finite element modeling in 
Comsol Multiphysics, the Seebeck voltage that is gener- 
ated due to the heating of the ferromagnet. In this model 
we impose the constant heat fiux, given by Eq. [4j through 
the top layer of the ferromagnet and solve the thermo- 
electric model [ini HZ] . Both, the heat flux and the calcu- 
lated Seebeck voltage are plotted in Fig. [3j: (solid lines) 
for multiple frequencies. Peak heights ranging from 98 
till 197 nV are calculated. For comparison, the observed 
peak height of Fig. |3^ is replotted in Fig. ^ as black 
dots. 

For a fixed rf field strength, the Seebeck voltage should 
increase monotonically with frequency due to an increas- 
ing dissipation (see Eq. |4|. However, because of the ex- 
perimental variation in rf field strength for different fre- 
quencies, a specific relation between the calculated See- 
beck voltage and the frequency is found (Fig. [3j:) . The 
experimental data is in agreement with the calculations 
within a factor of two and follows partially the same 
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tion dynamics and caloritronics. 
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FIG. 4. (color online) a) AMR voltage vs. magnetic field. 
Peaks and dips are observed for the resonance fields of the 
two driving frequencies (5 GHz apart). The different traces 
are offset for clarity reasons. b) The magnitude of the 
rf field is extracted from the peak height of the resonance 
{V = 7dcAi?sin^ ^c). c) Anisotropic magnetoresistance mea- 
surement with I and B perpendicularly aligned. 



trend. This discrepancy is attributed to small sample 
to sample variations in the rf field strength. Since the 
AMR measurements are performed on a separate device, 
the fields can differ for the thermoelectric device. More- 
over, small shifts in the contact area of the thermocouple 
can lead to changes in the heat transport and hence, dif- 
ferent thermo voltages. 

Furthermore, circulating rf currents combined with an 
oscillating magnetoresistance at the same frequency can 
cause dc voltages via a rectifying effect and mimic the 
observed thermal behavior in our devices |15| . We have 
excluded these effects by using a similar device with an 
Au-Au thermocouple such that the Seebeck effect van- 
ishes. For this device we observed a flat background 
voltage without peaks and dips. We note that thermal 
voltages can be of importance in other device geometries 
where ferromagnets are electrically connected to nonmag- 
netic metals. For example, detection of interface voltages 
that arise due to spin-pumping |19[ 120) cannot be easily 
distinguished from generated Seebeck voltages. Thor- 
ough temperature or material dependent measurements 
might offer a solution to discriminate between both mech- 
anisms. 

In conclusion, we have demonstrated a new thermo- 
electric detection technique for ferromagnetic resonance. 
The observed resonance peaks are in good agreement 
with the LLG equation and thermoelectric finite element 
modeling. Additionally, this technique can be applied on 
the nanoscale and is not limited to conductive ferromag- 
netic media. Thermal detection offers a valid, alterna- 
tive method for studying the dissipation, i.e. the Gilbert 
damping term, in nanoscale ferromagnetic islands. We 
hope that these results stimulate research for new physi- 
cal effects that arises from coupling between magnctiza- 



* 'f.l.bakker@rug.nl 
[1] K. Uchida, S. Takahashi, K. Harii, J. leda, W. Koshibae, 
K. Ando, S. Maekawa, and E. Saitoh, [Nature 455, 778| 
(20 08) 

[2J C. M. Jaworski, J. Yang, S. Mack, D. D. Awschalom, 
J. P. Heremans, and R. C. Myers, Nature Materials 9,| 
89 8 (2010) 

[3J K. Uchida, J. Xiao, H. Adachi, J. Ohe, S. Takahashi, 
J. leda, T. Ota, Y. Kajiwara, H. Umezawa, H. Kawai, 
G. E. W. Bauer, S. Maekawa, and E. Saitoh, Nature 



[12 

[13: 

[14 
[15- 



[17 
[18 
[19 
[20 

[ 



and 



Materials 9, 894 (2010) 

J. Xiao, G. E. W. Bauer, K. C. Uchida, E. Saitoh 
S. Maekawa, Phys. Rev. B 81, 214418 (2010) 

H. Adachi, J.-i. Ohe, S. Takahashi, and S. Maekawa, 
|Phys. Rev. B 83, 094410 (2011), 

S. I. Kiselev, J. C. Sankey, I. N. Krivorotov, N. C. Emley, 
R. J. Schoelkopf, R. A. Buhrman, and D. C. Ralph, 
Nature 425, 380 (2003) 

D. Houssameddine, U. Ebels, B. Delaet, B. Rodmacq, 

I. Firastrau, F. Ponthenier, M. Brunet, C. Thirion, J.- 
P. Michel, L. Prejbeanu-B uda, M.-C . Cyrille, O. Redon, 
and B. Dieny, Nature materials 6, 441 (2007) 

A. Brataas, Y. Tserkovnyak, G. E. W. Bauer, and 

B. I. Halperin, Physical Review B 66, 060404 (2002), 

Y. Tserkovnyak, A. Brataas, and G. E. W. Bauer, |Phys.| 
Rev. B 66, 224403 (2002) 



Y. Tser kovnyak, A. Brataas, and G. E. W. Bauer, |Phys. 
Rev. Lett. 88, 117601 (2002), ' 



A. Brataas, Y. Tserkovnyak, and G. E. W. Bauer,, Phys. 



Rev. Lett. 101, 037207 (2008), 

R. Meckenstock, Rev. Sci. lustrum. 79, 041101 (2008)| 
F. Sakran, A. Copty, M. Golosovsky, D. Davidov, and 
P. Monod, Applied Physics Letters 84, 4499 (2004) 
J. Coey, Magnetism and magnetic materials (Cambridge 
University Press, 2010). 

M. V. Costache, S. M. Wat ts, M. Slad kov, C. H. van der 
Wal, and B. J. van Wees, [Applied Physics Letters 89, | 
232115 (2006) 

F. L. Bakker, A. Slachter, J. - P. Adam, and B. J . van 
Wees, Physical Review Letters 105, 136601 (2010)| 
A. Slachter, F. L. Bakker, and B. J. van Wees, ArXiv 
e-prints (2011), arXiv: 1107.3290 [cond-mat .mes-hallj 
C. Kittel, Introductton to Solid State Physics -7th ed. 
(John Wiley & Sons, Inc., New York, 1995). 
X. Wang, G. E. W. Bauer, B. J. van Wees, A. Brataas, 
and Y. Tserkovnyak, [Phys. Rev. Lett. 97 , 216602 (2006)] 
M. V. Costache, M. Sladkov, S. M. Watts, C. H. van der 
Wal, and B. J. van Wees, Physical Review Letters 97^ 
^6603 (2006) 



